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The contribution deals with the experimental and numerical investigation of compressible flow through the 


tip-section turbine blade cascade with the blade 54” long. Experimental investigations by means of optical (inter- 


ferometry and schlieren method) and pneumatic measurements provide more information about the behaviour and 


nature of basic phenomena occurring in the profile cascade flow field. 


The numerical simulation was carried out by means of the EARSM turbulence model according to Hellsten [5] 


completed by the bypass transition model with the algebraic equation for the intermittency coefficient proposed 


by Straka and Příhoda [6] and implemented into the in-house numerical code. The investigation was focused par- 


ticularly on the effect of shock waves on the shear layer development including the laminar/turbulent transition. 


Interactions of shock waves with shear layers on both sides of the blade result usually in the transition in attached and/ 


or separated flow and so to the considerable impact to the flow structure and energy losses in the blade cascade. 
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Introduction 


Recent advances in steam turbine efficiency have lead 
to steam turbines with relatively long rotor blades of the 
last stages in low pressure cylinder. Considerable length 
of these blades results from pursuit of larger exit annulus 
area since larger exit annulus area leads to lower exit 
losses due to decreased exit velocity. Also, increasing the 
exit area leads to lower number of turbine rows and thus 
to lower production costs. However, aerodynamic design 
of efficient long blades is a challenging task since there is 
not much space for the design which often results in 
blades with tips operating under supersonic regimes [1], 
[2], [3], [4]. Therefore, both experimental and numerical 
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investigations of flow past blade cascades representing 
such supersonic tip sections are essential for successful 
design. 

Recent investigations showed that among other phe- 
nomena attention has to be paid to the interaction of 
shock waves with either suction or pressure side bound- 
ary layer. Namely interaction of front shock wave inner 
branch with either pressure side boundary layer or near 
wake structure significantly contributes to energy losses 
[1]. In order to understand and predict behaviour of these 
phenomena, detailed experiments and advanced numeri- 
cal simulations has to be performed. Such numerical 
simulations have to take into account many factors, espe- 
cially the effect of the pressure gradient, free stream tur- 
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bulence, and possibly wall roughness. The mathematical 
model should not be based on the turbulence model only, 
but should be completed by adequate models of the 
laminar/turbulent transition and turbulent heat transfer. 
Moreover in transonic and supersonic flows, the interac- 
tion of shock waves with shear layers on blade surface 
may cause the transition to turbulence and possibly the 
flow separation. 

The contribution deals with the numerical simulation 
of 2D compressible flow through supersonic turbine 
blade cascade with profiles corresponding to the tip- 
section of the low-pressure part of the steam turbine. 
Simulations carried out by the in-house numerical code 
based on the finite-volume solver for the RANS equa- 
tions closed by the explicit algebraic model Reynolds 
stresses completed by the modified algebraic transition 
model. Simulations for the supersonic inlet conditions 
were focused especially on the effect of shock waves on 
the development of shear layers on blades including the 
laminar/turbulent transition in attached and/or separated 
flows. Behaviour of the investigated phenomena is also 
studied and documented by classical optical methods 
such as interferometry and schlieren technique. 


Investigated Blade Cascade 


The tip section of the last rotor wheel, comprising of 
titanium rotor blades 54" long, is characterized by very 
large stagger angle of the tip section profiles. The con- 
vergent - divergent tip profile cascade is designed for a 
very small flow turning angle. The fan of these long 
blade airfoils is highly spaced on the tip diameter (the tip 
section will rotate on a diameter of 4700 mm), so the 
overlapping of the adjacent profiles is very small (Fig. 1). 

The profile applied to the tip section (Fig.1) is rela- 
tively very thin, and its central line is unusually cam- 
bered, so that shapes of the suction and pressure sides 
have converse curvatures: the nature of the shape of the 
suction side is therefore concave and, conversely, the 
shape of the pressure side is convex. Characteristic di- 
mensions of the blade cascade are shown in Table 1. 


Fig. 1 The convergent-divergent tip profile cascade 
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Table 1 Characteristic dimensions of the blade cascade 


Pitch/Chord s/b 0.951 
Max. Thickness/Chord T/b 0.025 
Throat Opening/Chord olb 0.169 
Stagger Angle y 78.37° 
Trailing edge Thickness/chord Trg/b 0.0033 
Design Inlet Flow Angle Oides 81.25° 


The profile cascade consists of 8 prismatic blades only. 
The profile chord of the model was chosen b= 0.15 m, 
i.e. the reduction scale is approximately 1:2, considering 
the real blade dimensions. The blade height is /= 0.16 m, 
so the aspect ratio is AR = 1.07 only. The model dimen- 
sions result from many requirements and restrictions (e.g.: 
corresponding Reynolds number values, solidity of the 
blades, dimensions of the test section, etc.). 


Experiment 


The experiments were carried out in the Aerodynamic 
laboratory at Novy Knin. The high-speed intermittent 
blow-down wind tunnel for 2D cascade testing was used, 
and the flow medium was dry air. A special test section 
was used for an investigation of the flow in a prismatic 
profile cascade with a small flow turning angle and su- 
personic inlet velocity. This test section has been de- 
scribed in [5]. 

The high-speed wind tunnel for 2D cascade measure- 
ments (see Fig. 2) is equipped with an adjustable super- 
sonic inlet nozzle. The nozzle has parallel side walls and 
the shape of the upper and lower wall can be continu- 
ously changed, so that inlet Mach number up to M; = 2.0 
can be reached. The air flows from the test section to the 
settling chamber. The backpressure is controlled by an 
adjustable control nozzle, which is situated downstream 
the settling chamber. The scheme of the test section is 
depicted in Fig. 3. Both the upper wall and the lower wall 
of the inlet channel are equipped with perforated inserts. 
These inserts are connected to the low pressure parts of 
the wind tunnel, and the suction is controlled by slide 
valves. The suction mode of the lower perforated wall is 
used for inlet velocities around a value of M, = 1.0, to 


1. silicagel dryer, 2. filters, 3. entrance nozzle, 4. inlet nozzle, 
5. transient insert, 6. rotatable test section, 7. settling chamber, 
8. control nozzle, 9. quick-acting valve, 10. diffuser, 11 main 
duct 


Fig. 2 Scheme of the suction type high-speed wind tunnel 
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Fig. 3 Test section with profile cascade arrangement for optical measurements 


improve the periodicity of the cascade inlet flow field. 
The upper perforated wall (with suction or with ventila- 
tion only) effectively prevents reflection of the outer 
branch of the inlet shock waves back into the studied 
flow field in the whole range of supersonic inlet Mach 
numbers, i.e. for 1<M,<2. The adjustable perforated 
tailboard is usually placed next to the trailing edge of the 
lateral profile in order to improve the periodicity of the 
exit flow field. In this specific configuration, the perfo- 
rated tailboard is situated downstream the blade next to 
the lateral blade. The optimal value for the angle between 
the tailboard and the plane of the trailing edges of the 
cascade was found experimentally. The tailboard not only 
prevents exit shock waves from reflecting back to the 
exit flow field, it also prevents expansion taking place at 
the lower wall block edge from influencing the exit flow 
field [5]. 

The inlet flow field parameters are measured by a 
number of static pressure taps on the sidewalls and by a 
Prandtl probe, which was placed upstream the cascade 
inlet. The relatively remote position of the probe is aimed 
at reducing the probe disturbances in the proximity of the 
leading edges of the cascade. The representative value of 
the inlet Mach number is evaluated from the average 
value of the static pressures, which is obtained from six 
pressure taps before the cascade (see Fig. 3) and the total 
pressure at the inlet (measured by the Prandtl probe), 
respectively. The isentropic exit Mach number was de- 
termined according to the static pressure, which was mea- 
sured in the settling chamber downstream the cascade. 

Optical measurements were carried out (interferome- 
try — infinite fringe method, schlieren method in the Toe- 
pler configuration). The field of vision for the optical 
measurements covered an important part of the cascade 
middle channel only (see Fig. 3). The dimension of this 
optical field is limited by the diameter of the glass win- 
dows for interferometry (180 mm). 

Experimental investigations were made at regime 
M,= 1.45 and Mis = 2.0. Corresponding Reynolds num- 
ber related to the profile chord b and to the isentropic exit 
Mach number M>;, in the wind tunnel (dry air) was dur- 
ing the experiment Re,,, = 1.92 x 10°. This value is of the 
same order as that of real operating turbine where the 
value is Re = 1.1 x 10°. The Axial Velocity Density Ratio 


(AVDR), which to some extent characterises the extent of 
3D phenomena, was monitored during the pneumatic 
measurements; the AVDR value for nominal condition 
was 1.059. 

Level of free stream turbulence is estimated from ear- 
lier measurements with transonic cascade to 1+2%. 


Mathematical model 

The mathematical model of compressible flow is 
based on the conditionally-averaged Navier-Stokes equa- 
tions closed by the EARSM turbulence model by Hell- 
sten [6] with the algebraic bypass transition model by 
Straka and Příhoda [7]. The EARSM model is modified 
into the form corresponding to models with the turbulent 
viscosity and so the turbulent stress is given by the relation 


7, ou, 7 
ty = Mh ou; , j S 
Ox; Ox; 3 ` OX, 


(1) 
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where the extra-anisotropy tensor a is given by 


relation 
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+ By (BSa Ann = QQ Sin ) 
The coefficients £; depend on invariants of strain-rate 
and vorticity tensors in the non-dimensional form. 
Transport equations for the turbulent energy k and the 


specific dissipation rate œ are given by the SST model 
according to Menter [8] in the form 
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Turbulent viscosity is given by the relation 


L = C, Pkr, (5) 
with the turbulent time scale 
T, = Max : xc, Y (6) 
po B ok 


where the Kolmogorov viscous time scale is used near 
the wall. The variable coefficient C, is obtained from the 
equation 
C,=-(B, +g Be)/2 (7) 
For the reduction of the undesirable over-production 
of the turbulent energy in the stagnation region, the pro- 
duction term in the turbulent energy equation was modi- 
fied into form 
P, = ve - 


i 
y ex; 


Jra y) uR? (8) 
where 7, is the Reynolds stress tensor and © is the 
absolute value of the vorticity tensor. The parameter y is 


given by the relation 
aha 
C, = 0.05 + 0.2 (9) 


The parameter wy is equal to 1 nearly in all the compu- 
tational domain except the vicinity of the stagnation point 
where it decreases to zero. The EARSM model is de- 
scribed in detail by Hellsten [6]. 

For the prediction of the transitional flows, the pro- 
duction and destruction terms in the k-equation are mul- 
tiplied by the intermittency coefficient y. Similarly, the 
effective viscosity is given by ef= 4+ yu, in the transi- 
tion region. The transition model is based on the concept 
of different values of the intermittency coefficient in the 
boundary layer y and in the free stream y,. The intermit- 
tency coefficient in the boundary layer 7; is expressed by 
the relation 


Sao 


y =1-—max | min (C7, 


7; =1-exp| -ño (Re, —Re,, ie (10) 


The transition onset is given according to Straka and 
Příhoda [7] by the empirical correlation for the momen- 
tum Reynolds number Reg = f (Tu, A) where Tu (%) is 
the free-stream turbulence level and A, is the pres- 
sure-gradient parameter. The length of the transition re- 
gion is expressed using the parameter N= 
where ñ is the spot generation rate and ø is the spot 
propagation rate introduced by Narasimha [9]. The effect 
of the free-stream turbulence and the pressure gradient on 
the parameter N is correlated by an empirical relation 
proposed for the attached flow by Solomon et al. [10]. 
The onset of transition in separated flow is given by the 
correlation proposed by Mayle [11] in the form 
Re, =f (Reg, Rex) where Reg is the momentum Rey- 
nolds number at the separation and Re,, is the Reynolds 


fio Reg 
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number related to the distance of the separation from the 
leading edge. The length of the transition region is ex- 
pressed according to Mayle [7] by the relation 
„228x10 
no = ——— 
Reg; 

To avoid the application of local variables, the maxi- 
mum of the vorticity Reynolds number Remax is used 
instead of the momentum Reynolds number Reg. This 
link is expressed by the relation Reg = Re max’ C where the 
parameter C depends on the pressure gradient. The used 
version of the algebraic transition model is described in 
detail by Fiirst et al. [12]. 

The algebraic model was implemented into the in- 
house numerical code. The code is based on the finite- 
volume method of the cell-centered type with the Osher's- 
Solomon's approximation of the Riemann solver and a 
two-dimensional linear reconstruction with the Van Al- 
bada's limiter. The governing equations are discretized 
using a multi-block quadrilateral structured grid with a 
block overlapping implementation. 


(11) 


Calculation setup 


Due to supersonic inlet conditions, the flow through 
the blade cascade is influenced by parasitic shock waves 
arising by the reflection from the computational domain 
boundaries. The application of the quadrilateral block- 
structured computational grids with their overlapping 
leads moreover to the reflection of shock waves on the 
block boundaries. Therefore the computational domain 
was adequately prolonged before and behind the blade 
cascade to weaken these reflections. The detail of the 
computational grid of the "chimera" type near the leading 
and trailing edges is given in Fig. 4. 


Fig. 4 Detail of the computational domain "chimera” 


The constant values of the total pressure po), total 
temperature Tọ; and incidence angle œ; were prescribed 
as inlet conditions. The static pressure was extrapolated 
from the computational domain. The static pressure given 
by the outlet isentropic Mach number was prescribed as 
the outlet condition. Inlet free-stream turbulence was 
chosen Tu = 5%. The ratio of viscosities 44/12 was chosen 
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according to the length of inlet part. The free-stream tur- 
bulence in the distance of one spacing before the cascade 
was Tu = 1.2%. 


Results 


The flow field in the interblade channel is supersonic, 
with a small subsonic region around the leading edge. 
The sonic throat is missing. In the relatively small area of 
the interblade channel, very intense supersonic simple 
expansion takes place on the suction surface near the 
leading edge (Fig. 5). 

Figure 6 compares the distributions of the isentropic 
Mach number along the suction side and the pressure 
side of the profile. The first distribution is obtained from 
experiment, while the second distribution is obtained 


from CFD simulation (commercial code) described in [1]. 


The calculated flow field in this case is completely tur- 
bulent (Fig. 7). This results in a less complicated interac- 
tion of the boundary layer with the internal branch of the 
inlet shock wave on the pressure side of the profile. By 
contrast, the test result (Fig. 5) shows that the interaction 
is more complex and local separation takes place. This is 
because the boundary layer on the pressure side is proba- 
bly laminar up to the interaction area. The impact of this 
complex interaction is that it suppresses the expansion 
area on the suction side of the neighbouring profile. The 
expected additional supersonic expansion upstream the 
point where internal branch of exit shock wave interacts 
with suction side boundary layer is completely cancelled 
(in the range 0.4 < x/bx < 0.6 for CFD and 0.2 < x/bx < 


Fig. 5 Interferogram showing subsonic region and expansion 
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Fig. 6 Distribution of isentropic Mach number along profile 
obtained by commercial code and experiment 
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Fig.7 Distribution of Mach number - fully turbulent CFD 


0.4 for experiment as can be seen from Fig. 6 ). There- 
fore, downstream this interaction, the values of Mach nu- 
mber up to the trailing edge are lower than in the case of 
turbulent interaction on the pressure side. This is also the 
reason why the interaction of exit shock with the suction 
side boundary layer is located further downstream in case 
of CFD (Fig. 6). This fact has an essential negative in- 
fluence on the aerodynamic loading of the studied profile. 


Current numerical results 


The distribution of the pressure coefficient along the 
blade obtained by the EARSM model with the algebraic 
transition model is shown in Fig. 8. The coordinate s is 
measured along the blade surface from the leading edge. 
The effect of the interaction of shock waves with bound- 
ary layer on the suction side is well visible at s/Smax % 0.2; 
0.35 and 0.5 approximately. This distribution does not fit 
that obtained from experiment perfectly, but it is clear 
that location of the exit shock/suction side boundary 
layer interaction is at S/Smax 0.5 unlike in case of fully 
turbulent calculation (Fig. 6, [1]). 

The field of Mach number isolines in the blade cas- 
cade is shown in Fig. 9. Because of the stagger angle and 
inlet flow conditions, the concave side of the blade is the 
suction type while the convex side is the pressure type. It 
can be seen in Figures 9, 10 and 11, that the reflected 
oblique shock wave impinges on the suction side where 
the interaction with the boundary layer causes the bypass 
transition in attached flow in the distance s/Smax ~ 0.15 
from the leading edge. On the contrary, the impact of the 
front shock wave from the adjacent blade leads on the 
pressure side to the transition in the distance s/s, = 0.88 
from the leading edge. 

The distribution of the Reynolds number Reg related to 
the momentum thickness of the boundary layer and of the 
skin friction coefficient Cp related to inlet flow parame- 
ters is shown in Fig. 10. 


201706.00726v1 


chinaXiv 


ChinaXivA (ERAT 


Martin Luxa etal. Investigation of the Compressible Flow through the Tip-Section Turbine Blade Cascade with Supersonic Inlet 143 


1.0 


CFD EARSM model 
+% experiment 


0.5 f 
o : pressure side 
s 
0.0 : . 
se side 
t Je * 
-0.5 
0.0 


Fig. 8 Distribution of the pressure coefficient along the blade 
surface obtained by inhouse code and experiment 
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Fig. 9 Mach number isolines in the detail of the computa- 
tional domain 


The bypass transition on the suction side is evident 
from the distribution of Cz. The interaction of the bound- 
ary layer with inner branch of the exit shock wave can be 
shown on the distribution of Reg and especially Cp in the 
distance S/Smax % 0.5 from the leading edge. Following the 
smooth development of the laminar boundary layer, the 
separated-flow transition with a short separation region 
and following reattachment occurs on the pressure side 
before the trailing edge due to the interaction with the 
front shock wave from the adjacent blade. 

Comparison of the flow field by CFD and schlieren 
technique in the interaction region is shown in Fig. 11. 
Both CFD and experiment show structure typical for 
shock/boundary layer interaction with local separation 
region. Impinging front shock causes boundary layer to 


separate with separation point upstream of the interaction. 


Effective change of the flown surface at the separation 
point and reattachment point results in origin of two left 
running shockwaves. In case of experiment, the separa- 
tion region is noticeably longer and the point of front 
shock impingement is closer to the trailing edge. 

The field of turbulent energy in the prolonged compu- 
tational domain is shown in Fig. 12. The flow along the 
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Fig. 10 Distribution of boundary layer parameters along the 
blade surface 


Fig. 11 Comparison of flow in the region of interaction of the 
front shock and pressure side boundary layer - 
schlieren pictures 
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Fig. 12 Field of turbulent energy in the detail of prolonged 
computational domain 
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whole blade is attached with the exception of the short 
separated bubble before the trailing edge of the pressure 
side. This is in accordance with the turbulent energy field, 
where the growth of turbulent energy is more apparent on 
the suction side with the earlier transition. The substan- 
tially higher increase of turbulent energy occurs only in 
the wake behind blades. 

The distribution of the relative total pressure p29/p j9 in 
the wake in the distance of 25 mm behind the cascade 
exit plane is shown in Fig. 13. Besides the prediction by 
means of the EARSM model with algebraic transition 
model, the ~Re transition model of Langtry and Menter 
[13] implemented in the commercial code ANSYS Fluent 
was used for the comparison. The distribution of total 
pressure corresponds to the mean velocity profiles in the 
wake behind blades and is a measure of energy losses in 
the blade cascade. Besides the actual wake, the effect of 
the outer branch of the exit shock wave is shown on the 
distribution of the total pressure. 
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Fig. 13 Distribution of the relative total pressure behind the 
blade cascade 


Conclusions 


The paper shows that problem of shock wave bound- 
ary layer interaction plays important role in case of recent 
designs of supersonic tip sections for long last stage rotor 
blades. It has been shown that in case of tip section blade 
cascades, very short region before trailing edge has great 
influence on the flow within the interblade channel. More- 
over, in case of supersonic tip sections, the flow in this 
region is influenced by front shock wave. From this point 
of view it is crucial whether the boundary layer on the 
pressure side can withstand interaction with this front 
shockwave. Therefore, to predict flow through such blade 
cascades correctly, it is necessary to employ transitional 
turbulence model when RANS is to be used. Results 
published in the paper show that the EARSM model per- 
forms reasonably well. 
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